Swapping and entangling hyperfine coupled nuclear spin baths 
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We numerically study the hyperfine induced nuclear spin dynamics in a system of two coupled 
quantum dots in zero magnetic field. Each of the electron spins is considered to interact with an 
individual bath of nuclear spins via homogeneous coupling constants (all coupling coefficients being 
equal). In order to lower the dimension of the problem, the two baths are approximated by two 
single long spins. We demonstrate that the hyperfine interaction enables to utilize the nuclear baths 
for quantum information purposes. In particular, we show that it is possible to swap the nuclear 
ensembles on time scales of seconds and indicate that it might even be possible to fully entangle 
them. As a key result, it turns out that the larger the baths are, the more useful they become as 
a resource of quantum information. Interestingly, the nuclear spin dynamics strongly benefits from 
combining two quantum dots of different geometry to a double dot set up. 
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Introduction-Electron spins confined in semiconduc- 
tor quantum dots with an s-type conduction band, like 
for example GaAs quantum dots, experience decoherence 
through the spin-orbit interaction, and by the hyperfine 
interaction with surrounding nuclear spins. With respect 
to possible future solid state quantum computation sys- 
tems utilizing the electron spin as the qubipEl, these in- 
teractions act as a source of decoherence. Due to the spa- 
tial confinement of the electron spin in a quantum dot, 
the relaxation time T\ induced by the spin-orbit interac- 
tion is enhanced for low temperatures^. As the dephas- 
ing time Ti due to the spin orbit interaction turns out 
to be as long as the T\ time under realistic conditions^, 
the major source of decoherence in semiconduct or q uan- 
tum dots results from the hyperfine interactiorPSSl. For 
related reviews the reader is referred to RefsP^^. Simi- 
lar situations arise in carbon nanotube quantum dots^, 
phosphorus donors in silicon^ and nitrogen vacancies in 
diamoncPEl. 

Apart from this detrimental effect of the hyper- 
fine interaction, it provides a way to efficiently access 
the nuclear spins by e.g. external degrees of free- 
dom. This for example enables to bu ilt up an in- 
terface between light and nuclear spin d 21 l 22 J, to po- 
larize n uclear spin baths^H^S, to set up long-lived 
quantu m 126 ! 27 ! and classical memory devices or to gen- 
erate entanglementPSl. 

In both of the aforementioned contexts it is of key im- 
portance to understand the hyperfine induced spin dy- 
namics. Here one has to distinguish between the case of 
a strong and the case of a weak magnetic field applied to 
the electron spins. In the first limit, the "flip-flop" terms 
between the electron and the nuclear spins occurring in 
the Hamiltonian are strongly suppressed. This allows to 
treat them perturbatively or to even completel y negl ect 
them, which strongly simplifies the calculation :] 6 ! 7 ! 30 * ^. 
In the absence of such an external magnetic field, how- 
ever, many approximative techniques break down, and 
one has to resort to exact methods. As explained in 
RefP^, in order to gain exact results, strong restrictions 



on the initial state^^, the s ize of the systeirf^^ or the 
hyperfine coupling constants^SSSSS have to be made. 

In the present paper we combine the second and the 
third approach and focus on the advantages of the hy- 
perfine interaction. To this end we consider a model of 
two exchange coupled electron spins each of which is in- 
teracting with an individual bath of nuclear spins. This 
corresponds to the situation of spatially well-separated 
quantum dots. Assuming the baths to be strongly po- 
larized in opposite directions initially, we investigate, by 
means of exact numerical diagonalization, to what ex- 
tent it is possible to swap and entangle them. Usually 
exact numerical diagonalizations are restricted to rather 
small system sizes. In order to go beyond these limits, 
we reduce the dimension of the problem by approximat- 
ing the two baths by two single long spins. The spectral 
properties of the model described above have recently 
been studied in RefP^, where it has been shown that the 
spectrum of the Hamiltonian exhibits systematically de- 
generate multiplets. Motivated by these findings, below 
we will distinguish between an inversion symmetric sys- 
tem, showing the mentioned degeneracies, and a system 
with broken inversion symmetry where such degeneracies 
are absent. 

The work presented in this paper complements the re- 
sults of RefpS, where we analytically studied the homo- 
geneous coupling case for two electron spins coupled to 
a common bath of nuclear spins. Choosing the hyperfine 
coupling constants to be equal to each other of course has 
to be regarded as a rough approximation. However, it has 
been shown that already such simple models can y ield 
concrete predictions and realistic resultd 2 ^! 22 ! 24 ! 29 ! 3 ^. 

Model and methods. -The Hamiltonian of two exchange 
coupled electron spins, each of which is interacting with 
an individual bath of nuclear spins via the hyperfine in- 
teraction reads 

Ni N 2 

H = S 1 -Y / +S 2 -J2 A% + J eJi ■ S 2 , (1) 
i=l i=l 

where Sj are the electron and Zy are the nuclear spins 
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the j-th electron spin interacts with. For simplicity we 
will consider Ni — N 2 =: N in what follows. The param- 
eter J ex denotes an exchange coupling between the two 
electron spins, which can experimentally be adjusted in 
a range of [— 10 -3 , 10 _3 ]eV, and A], Aj are the hypcrfinc 
coupling constants. In a realistic quantum dot these are 
proportional to the square modulus of the electronic wave 
function at the sites of the surrounding nuclear spins. As 
a typical example, in GaAs quantum dots the overall cou- 
pling strength of the j-th electron spin A^ := $Z i=1 -M ^ s 
of the order of [l(T 4 , lCT 5 ]eV. 

Due to the spatial variation of the electronic wave 
function, the hyperfine couplings are clearly inhomoge- 
neous. However, in the following we consider them to be 
equal to each other, meaning that A\ = A> /N. Then 
the Hamiltonian (TlJ) conserves apart from the total spin 



J = Si + S 2 + h + h, where I 3 
squares of the total bath spins IJ 



E 



1 Iij , also the 



H,J 



H,I* 



= 0. 



(2) 



The first symmetry will be helpful for the e xact n umerical 
diagonalizations of the Hamiltonian matrbpEID through 
which we will obtain the dynamics in the following. We 
compute the time dependent density matrix by decom- 
posing the initial state into energy eigenstates and ap- 
plying the time evolution operator. Tracing out the elec- 
tron degrees of freedom then yields the reduced density 
matrix of the nuclear baths p n (t) from which we can cal- 
culate the time evolution of all observables. For details 
the reader is referred to Ref!^. 

In the following we will approximate each of the two 
baths Ij by a single long spin. Let us briefly discuss 
to which physical situation this corresponds. A general 
state of a bath is a superposition of states from different 
multiplets 



E 



(3) 



where the quantum numbers associated with a certain 
Clebsch-Gordan decomposition of the bath have been 
omitted. The number of multiplets contributing to the 
sum in ([3| decreases with increasing bath polarization. 
For very high polarizations we can therefore approximate 
the state Q by \fij) — \I, rrij) with mj ±1. Due to the 
commutation relations ^ all dynamics is then captured 
by the following simple Hamiltonian, to which we refer 
to as the long spin approximation Hamiltonian: 



LSA 
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Si-h 



21 



S 2 -h + Jex Si ■ S 2 (4) 



The form of the couplings A 5 /2I results from the ob- 
servation that the N bath spins can couple to I = 
N/2, N/2-1, N/2-2, .... As we assume highly polarized 
baths, we consider the maximal value / = N/2. Solving 
for N then yields the coupling constants in Q . For later 
convenience we define A := A 1 + A 2 . 
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FIG. 1: Nuclear spin dynamics of the Gaudin model for 
inhomogeneous (upper panels) and homogeneous (lower pan- 
els) hype rfine c ouplings. The bath state is initially randomly 
correlate d 11 * 35 ! and the electron spin is pointing upwards. The 
left column shows the case of N = 30 and a bath polarization 
of 0.93, while in the right column we have N = 12 and a po- 
larization of 0.33. Even in the latter case the bath dynamics 
for inhomogeneous and homogeneous couplings are still quite 
similar to each other. 



High nuclear polarizations of up to 80% have been ex- 
perimentally demonstrated in RefsP^H. In particular, 
Refp3 also discussed the possibility of polarizing two nu- 
clear ensembles in different directions. However, a ques- 
tion concerning the LSA arises from assuming the cou- 
plings to be homogeneous: As demonstrated in Rcf. 34 , 
this approximation is a good one for short time scales, 
whereas for longer times artifacts occur. As the nuclear 
dynamics are slow, it has to be questioned to what extent 
homogeneous couplings are adequate in order to evaluate 
nuclear spin dynamics. Therefore we numerically inves- 
tigated the time evolution of the nuclear spins in the 
Gaudin model. The Gaudin model is the central spin 
model with a single central spins, as corresponding to 
one of the two first terms in the Hamiltonian (TJ). As a 
result from our numerics, we find that the influence of 
inhomogeneities is suppressed with increasing polariza- 
tion. This is illustrated Fig. [l] where we compare cases 
of high and low polarization. Even in the latter case the 
dynamics for both types of couplings are very close to 
each other. In the LSA we are considering two coupled 
Gaudin models. Consequently, the data presented in Fig. 
[TJdoes not give a strict proof for the adequacy of neglect- 
ing inhomogeneities in the couplings. However, it clearly 
provides qualitative evidence into this direction. 

Swapping nuclear spin polarizations. -We will now eval- 
uate the nuclear spin dynamics within the LSA and ex- 
plore the possibility of swapping oppositely polarized spin 
baths, (II) = —(Iij,), Our initial state \a) will be a 
simple product state between the electron and the nu- 
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FIG. 2: The nuclear baths are considered to consist of spins 
with length (1/2). The frame marks the LSA system with 
the smallest possible bath spin length / = 1/2. Here the 
dynamics of all four spins are highly coherent provided the 
values of all couplings are close to each other, motivating the 
condition «5J> (see text). 
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FIG. 3: Nuclear spin dynamics for J ex /(A/2I) = 3.5, A = 
and \a) = (1/VT3) (2|1W + 3|M)) W - 1,1). The upper 
panel shows data data I — 80 and M = where the expecta- 
tion values (Ij(t)) decay rapidly to zero. In the bottom panel 
we consider / = 1500 and M = 159 where an almost complete 
swap of the nuclear spins is observed. 



clear state \a) = |a e )|a„). Since both baths are spa- 
tially well-separated, the initial nuclear state is again a 
product state of the two long spins. Note that entan- 
glement within the baths can not be considered within 
the LSA. However, due to their high polarizations, cor- 
relations within the baths are of minor importance. In 
what follows, we will always work in subspaces of fixed 
J z =: M where only the z-component has a non-zero 
expectation value (7J(t)). Note that (again) due to the 
high bath polarizations, it is realistic to assume that the 
initial state has components exclusively in subspaces of 
fixed M. Moreover, we will assume the z-component of 
the total electron spin to be initially zero, i.e. the spins 
are antiparallel. Similar results as to be presented below 
are obtained for more general initial states of the electron 
spin system. 

More importantly, we will concentrate on exchange 
couplings being of the same order of magnitude as the 



hyperfine coupling strength, 

J ex /{A/2I) a 1, (5) 

meaning that the two electron spins are coupled as 
strongly to each other as they are coupled to the bath 
spins. This is motivated by the following observation. 
Let us consider ^ for the smallest possible value, I = 
1/2. As shown by elementary numerics, the dynamics 
of all four spins are, under the condition (J5|, highly co- 
herent, and the nuclear spin polarizations can nicely be 
swapped, i.e. at the end of the process the expectation 
values (Ij{t)) are, to a very good degree of accuracy, ex- 
changed as compared to the initial state. Let us now 
consider the two baths in the original model ([I]) to con- 
sist of spins with length (1/2), as already assumed for 
the derivation of the couplings in Q. As depicted in 
Fig. [2j the complete system can now be regarded as set 
of I = 1/2 models. Thus, from a heuristic point of view, 
the biggest chance to swap the full baths exists if all the 
subsystems are swapped. Hence, the exchange coupling 
has to be of the order of the coupling between the electron 
and the bath spins for any subsystem. For homogeneous 
couplings within the baths, this means that J ex ps A/N, 
which translates into ^ for the LSA, as explained above. 

However, at first sight, it does not seem to be possible 
to swap the initially antiparallel nuclear spins I > 1/2, 
even if the condition ^ is fulfilled. This is demon- 
strated in the upper panel of Fig. |3j where we con- 
sider I = 80, J ex /(A/2I) — 3.5 and a zero "detuning" 
A := A2 — A\ = 0. This corresponds to a situation 
in which the two quantum dots have the same geome- 
tries. We choose the comparatively generic electron spin 
state \a e ) = (l/\/l3) (2|-ft\|J,) +3|M)) (similar results oc- 
cur for other choices) and plot the dynamics for antiparal- 
lel nuclear spin configurations with the maximal possible 
z components \a n ) = As seen from the figure, 

the expectation values of the bath spins decrease quite 
rapidly to zero and are far away from being properly 
swapped. 

Surprisingly, this apparently negative result turns out 
to be related to the size of the baths in the following 
sense: Consider an initial state with the electron spins 
being in an arbitrary linear combination of Iff-ll), l-llft") 
and a nuclear state fulfilling, say, 1\ < and 

7| > with \If\ < |/f I, i.e. the "magnetization" M r := 
M t '(21 + 1) is nonzero, M r = (|Jf | - \I!\)/{2I + 1) ^ 0. 
Here we find that for magnetizations larger than a certain 
"critical" value M^r (slightly depending on the electron 
spin state) the expectation value {If{t = 0)) is, to a 
an excellent degree of accuracy, completely reversed. As 
a representative example, in the left panel of Fig. [4] we 
plot the quantity for I = 200 with, as before, 

\a e ) = (I/a/13) (2|1HJ.) + 3|^)) as initial electron spin 
state. We consider two different initial nuclear states 
\a n ) = \M — I, I), where the corresponding value of M r 
is in one case exactly at, in the other case lower than the 
critical Af^r at the given nuclear spin length. In the latter 
case the reversal of (I[(t)) is slightly incomplete. 



4 



-M=M r =0.11 
M r =0.06 

r 




t [10 hI/(7t)A] 



FIG. 4: Left panel: (/'(*)) for J ex /(A/2I) = 3.5, 
A = and / = 200. The initial state is \a) = 
(1/a/13) (2|fJJ.) + 3|M>) |M - 1,1} where two values of the 
magnetization M r — M/(2I + 1) are considered; M r = M^. = 
0.11 corresponds to the critical value at I — 200. The horizon- 
tal lines are guides to the eye indicating the value needed for 
a complete reversal of {11(f)}. Right panel: Af^r versus spin 
length I for the same initial state as in the left panel: The fit 
results in a power-law decrease M°(I) = e -°- 5 j-°- 33 Thus, 
for large enough spin baths antiparallel nuclear spin configu- 
rations can, to an excellent degree of accuracy, be swapped. 
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FIG. 5: Left panel: Reversal of (If(i)} for the same situa- 
tion as in Fig. [4] left panel, but A = -1.2A and M r = 0.09: 
Breaking inversion symmetry facilitates swapping the nuclear 
spin polarizations. Right panels: Swap time r at A = and 
Jex/(A/2I) = 3.5 as a function of / (upper panel), and at 
/ = 200 as a function of (J ex /(A/2I)) (lower panel). In the 



upper case we find a power law as r = e 



0.015 rl.19 



the lower case r = e 7 ' 51 [J ex /(A/2I)]~ L - m6 (hl/nA). 



(hl/nA), 



It is now a key observation that Mf(I) strongly de- 
creases with increasing spin length /. This is demon- 
strated in the right panel of Fig. [4] for spin lengths up to 
I = 600, where a clear power law scaling is found: 
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(6) 



Hence, for a large enough value of / the magnetization 
M^r will be so close to zero that, up to irrelevant correc- 
tions, antiparallel nuclear spin configurations can indeed 
be swapped. To give a quantitative example, for I = 10 6 
(as typical for the nuclear spin bath of GaAs quantum 
dots) the above power law leads to M° = 0.006 implying 
that baths with If = — 0.988/f can be swapped. Finally, 
in the bottom panel of Fig. |3]we plot (If(t)) and (Ii(t)} 
for I = 1500 (which is about the largest system size ac- 
cessible to our numerics) : Obviously we are very close to 
a full swap. 

Moreover, the performance of such a swap process can 
significantly be further improved by departing from the 
symmetric case A = A\ — A2 = 0, i.e. considering differ- 
ent geometries for the two quantum dots: The left panel 
of Fig.[5]shows (If(i)) for the same situation as in the left 
panel of Fig. |IJ but A = -1.2 A and M r = 0.09 (which 
is lower than M r c = 0.11 found before for A = 0). As 
seen, (If (t = 0)) is still fully reversed. Interestingly, this 
result turns out to be rather independent of the precise 
value of A ^ (including its sign), suggesting that the 
observed increase of "swap performance" goes back to 
some qualitative change in the dynamical properties. In 
fact, as shown in Ref!^, the spectrum of inversion sym- 
metric systems exhibit a macroscopically large subspace 
of energetically degenerate multiplets. Although the ini- 
tial states considered throughout this manuscript lie in 
energy quite far away from those degenerate levels, it is 
an interesting question to what extent both observations 
are related. 

Finally, in the right panels of Fig. [5] we analyze the 
duration t of the swap process as a function of the spin 
length / as well as the ratio J ex / (A/21) for again A = 0. 
In both cases we find power law dependencies leading for 
a realistic system size of I = 10 6 to a swap time of r of 
a few ten seconds. 

It is well-known that the nuclear bath is not static. 
The nuclear spins are interacting through e.g. dipolar 
and quadrupolar interactions^. These could be possible 
limitations to the phenomena described above. In order 
to circumvent the resulting problems, one would have to 
use additional techniques like e.g. refocusinj^l. 

Entangling the nuclear baths — In order to measure the 
entanglement between the long bath spins , we utilize the 
(logarithmic) negativity L defined by * 41 * 42 * 



£ = log 2 (||^| 



(7) 



where denotes the trace norm P||i = Tr(VA+A), 
and p\ is the partial transpose of p n with respect to the 
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first spin I\. Since 



\pnh 



(8) 



where Ef denote the eigenvalues smaller than zero, the 
negativity essentially measures to what extent the par- 
tial transpose fails to be positive, indicating non-classical 
correlations^. 

In order to evaluate the dynamics of the negativity, 
p x n has to be diagonalized in each time step considered. 
This is a numerical effort which restricts us to system 
sizes somewhat smaller than considered before. The left 
panels of Fig. [6] show the entanglement dynamics for two 
spin lengths / = 20, 80 at comparatively high polariza- 
tion M r = 0.8 and A = 0. The initial state is the same 
as used before, \a) = (1/-/13) (2\U) + 3|M» \M-I,I). 
In both cases the dynamics are rather similar to each 
other: A rapid increase of the negativity is followed by 
a more or less regular oscillation around a mean value 
which increases with the spin length /. In particular, the 
negativity never returns to zero. 

In order to quantify these observations we introduce 
a relative negativity L r = L/L max , where L max — 
log 2 (21 + 1) is an upper bound of L (cf. RefP^), and an- 
alyze the maximum L^P° of this quantity attained within 
a fixed interval [0, 200](hl /irA). The results are plot- 
ted in the right panel of Fig. [6j While the spin lengths 
achievable here are too small to allow for a quantitatively 
meaningful fit, the data still shows a significant growth 
with increasing I (suggesting, in fact, a power law). This 
observation implies that, similarly as for swapping nu- 
clear spin polarizations, also entangling spin baths bene- 
fits from large bath sizes. We note that this effect is not 
due to the simple growth of the reduced density matrix 
with increasing / since we are considering the relative 
negativity where such influences are scaled out. On the 
other hand, by the same argument, the maximal relative 
negativity should decrease with increasing magnetization 
at fixed /; an example for this behavior is shown in the 
left panel of Fig. [7] 

In the right panel of Fig. [7] we finally demonstrate 
the influence of a non-zero detuning for different spin 
lengths and magnetizations. Similarly to the results re- 
garding a nuclear swap, the entanglement is enhanced by 
a non-zero detuning with its precise value being again of 
minor importance. This supports the conjecture that the 
systematic degeneracy reported in Refpa has a clear dy- 
namical signature. Interestingly, breaking the inversion 
symmetry has stronger influence for higher magnetiza- 
tion. 

Conclusions. -hi summary we have studied the spin 
and entanglement dynamics of the nuclear baths in a 
double quantum dot. Each of the two electron spins was 
considered to interact with an individual bath via ho- 
mogeneous couplings. In order to lower the dimension 
of the problem, both baths have been approximated by 
long spins. We focused on the virtue of the hyperfine in- 
teraction and regarded the electron spins as an effective 
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FIG. 6: Left panels: Negativity L(t) for I = 20,80 and 
M r = 0.8, J ex /(A/2I) = 3.5, A = 0. Right panel: Maximal 
relative negativity L^ 00 in the time interval [0, 200](ft//7r^4) as 
a function of spin length / for otherwise identical parameters. 
We find a clearly increasing curve, indicating that for large 
enough sizes, the baths can be fully entangled. 
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FIG. 7: Left panel: Maximal relative negativity L^ 00 versus 
M r for I = 15, A = and the same initial state as in Fig. [6] 
Right panel: Lf )0 as a function of the detuning A for various 
parameters. The precise value of the detuning is again of no 
particular importance. 



coupling between the baths. We demonstrated that it is 
possible to swap them if their size is large enough, and 
provided strong indication that, under the same condi- 
tions, it might be even possible to fully entangle them. 
Surprisingly, it turns out to be advantageous to use dots 
of different geometry (enabling for A ^ 0) to built up 
the double quantum dot. 
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